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ITERATED BROWNIAN MOTION IN AN OPEN SET 

By R. Dante DeBlassie 

Texas A&M University 

Suppose a solid has a crack filled with a gas. If the crack reaches 
the surrounding medium, how long does it take the gas to diffuse out 
of the crack? Iterated Brownian motion serves as a model for diffusion 
in a crack. If r is the first exit time of iterated Brownian motion 
from the solid, then P{t > t) can be viewed as a measurement of the 
amount of contaminant left in the crack at time t. We determine the 
large time asymptotics of P(r > t) for both bounded and unbounded 
sets. We also discuss a strange connection between iterated Brownian 
motion and the parabolic operator ^A^ — 

1. Introduction. Suppose an infinite slab with finite thickness is cracked. 
The crack is fihed with a gas and the crack reaches the surrounding medium. 
How long does it take the gas to diffuse out of the crack? Burdzy and 
Khoshnevisan [14] show a reasonable model of diffusion in a crack is given 
by iterated Brownian motion. In essence, the crack is modeled as the path 
of two-sided Brownian motion. In analogy with ordinary Brownian motion 
and diffusion, if r is the first exit time of iterated Brownian motion from 
the solid, then P(r > t) provides a measure of the amount of contaminant 
left after time t. This is the object of our study. 

To define iterated Brownian motion Zt started at 2 S M, let and 
Yt be independent one-dimensional Brownian motions, all started at 0. Two- 
sided Brownian motion is defined to be 



Xt 




t >o, 

t <0. 



Then iterated Brownian motion started at z G M is 

Zt = z + X{Yt), t>0. 

In one requires X^ to be independent n-dimensional Brownian motions. 
Observe there is a choice here: ours or the requirement that the components 



Received June 2001; revised July 2003. 

AMS 2000 subject classifications. 60J65, 60K99. 

Key words and phrases. Iterated Brownian motion, exit time. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability^ 
2004, Vol. 14, No. 3, 1529-1558. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



R. D. DEBLASSIE 



of n-dimensional iterated Brownian motion be independent one-dimensional 
iterated Brownian motions. In fact, the latter process is the subject of the 
article [48] by Xiao. Our choice is motivated by a connection with — ^ 
and the interpretation of Burdzy and Khoshnevisan [14] as diffusion in a 
Brownian crack. 

Recently this process has been the subject of many articles. Limit the- 
orems are studied in [4, 17, 21, 29, 32, 46]; path properties and local time 
in [8, 11, 18, 30, 31, 48]; level sets in [13] and [28]. Burdzy [12] studied 
the variation of iterated Brownian motion. These results were extended by 
Khoshnevisan and Lewis [33] in their development of a stochastic integral 
for iterated Brownian motion. 

Although iterated Brownian motion is not a Markov process (the Chapman- 
Kolmogorov equation is not valid), a variant of it is known to have a con- 
nection with the parabolic operator — ^ [23]. 

The purpose of this article is two-fold: First, we establish a strange con- 
nection between iterated Brownian motion and a certain partial differential 
equation involving — ^. If this connection is analogous to the Brownian 
case, there might be a way to study the lifetime of iterated Brownian mo- 
tion in an open set by solving an initial-boundary value problem. We show 
this is not true in general. The second purpose of the article is to determine 
asymptotics of the distribution of the lifetime using another method. Now 
we write detailed statements. 

For typographical simplicity, write 

^"(*'")=(2^^"K"y- 

If / G C;,(M"), the space of bounded and continuous real- valued functions 
on R", set 

u{t,x)=EM{Zt)), 
where denotes expectation associated with Zq = x. Then 

n(t, x) = E(^f{x + X+(0 V Yt) + ^-(0 V {-Yt)))) 
(1.1) = I" i?(/(x + X+(0Vy)+X-(0V(-y))))pi(t,y)d2/ 



oo 
oo 



= 2/ / fiw)pniy,\w-x\)pi{t,y)dwdy. 

Jo JR" 

If we proceed formally and differentiate under the integral with respect to 
X, using that Apn{t,x) = 2^pn{t,x), we get 

poo r 

A'^u = 2 I I f{w)A,j.A,j:{pn{y,\w-x\))pi{t,y)dwdy 







oo 



JR"^ 



2Pn{yAw - X\) 



Pi{t,y)dwdy. 
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Next, change the order of integration, integrate by parts twice and then 
reverse the order of integration again to get 

A^u = 8 / f{w)pn{y,\w-x\)^-jpi{t,y)dwdy 
Jo JR" oy 

du 

Thus, there is hope a shght variation of Funaki's result carries through for 
our version of iterated Brownian motion. This is false. Differentiation under 
the integral is not possible and this leads to the following strange connection. 

Theorem 1.1. Let f £ C^(M"). Then the function u{t,x) = E^[f{Zt)] 
solves the Cauchy problem 

C7°°((0,oo) xM")nCfe([0,oo) xM"), 
(1.2) QA2-^^^,(i,x) = -i-i=A/(x) on(0,oo)xM-, 
u{0,x) = f{x). 

After submitting this article, an elegant variant of Theorem 1.1 due to 
Allouba and Zheng [2] appeared in print. See Theorem 0.1 as well as other 
interesting results in that paper. The theorem generalizes our Theorem 1.1 
in two ways: First, the Laplacian A can be replaced by the generator A 
of a continuous Markov process. Thus, iterated Brownian motion can be 
replaced by a more general iterated Markov process. Second, the function / 
need only be a bounded element of the domain of A with bounded Holder 
continuous second-order partial derivatives. 

Let Bt be re-dimensional Brownian motion and D C a reasonable open 
set. For td{B) = inf{t >0:Bt^ D}, it is well known that the function 

v{t,x) = EMBt)I{rDiB)>t)] 

solves the initial-boundary value problem 

1 d\ 
A- — h; = on(0,oo)xL>, 



.2 dt, 

f (0, x) = g{x), x £ D, 

v{t,x) = 0, x£dD. 

In light of Theorem 1.1, if D CI is an open set, there is hope the function 
u(t,x) = Px{td{Z) > t) solves 

^A2-^^n(t,x)=0 on(0,oo)xZ), 

u{0, x) = 1, X G D, 

u{t,x) = 0, x£dD. 
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Since the PDE is fourth order, an additional boundary condition is needed. 
For this choice of u, it turns out to be = 0, x € dD, where ^ is the 
inward normal derivative. Unfortunately, this fails. 

Theorem 1.2. Let D = (0, 1). Then for any a > 0, the function g{t,x) = 
Px{t(o,oo){Z) > t) does not satisfy 

d^g ^ dg_ 
^dx^ ~ dt' 

We use another method to study the lifetime of iterated Brownian motion 
in an open set D C R". For n-dimensional iterated Brownian motion Z, let 

TDiZ)=mf{t>0:Zt^D} 

be the lifetime of Z in D. Here and in the sequel we write f ~ g and f ^ g to 
mean for some positive Ci and C2, Ci < ^ < C2 and / < Cig, respectively. 

We also write f{t) ~ g{t) as t ^ 00 to mean — > 1 as t ^ 00. 

Let D C M" be an open cone, with vertex 0, such that 5"^^ H I? is regular 
for the Laplace-Beltrami operator L^n-i on S"'~^. We call D a general- 
ized cone. Then for some p{D) > 0, the exit time to{B) of n-dimensional 
Brownian motion from D satisfies 

Px{TD{B)>t)r.C{x)t~P^^^ ast^oo 

([5], Corollary 1). If n = 2 and the angle of D is ^ G (0,27r), then 

(1-3) ^'(^) = ^' 

see [15]. 

Theorem 1.3. Let D C M" he a generalized cone. Then as t —> 00 

P,{TD{Z)>t)^{t-Hnt, p{D) = l, 

[i-(p(D)+l)/2^ p(L>)>l. 

Since Brownian motion in a half-space has the same lifetime as Brownian 
motion in the positive reals, the next result is immediate from Theorem 1.5 
and formula (1.3). 

Corollary 1.4. If DQR has the form (—00, a) or {b,oo), then as 
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Theorem 1.5. Let D C R" be bounded and open with regular bound- 
ary. If Xf) is the principal eigenvalue of on D with Dirichlet boundary 
conditions, then 

logP^(Tz)(Z) >t) ~ -|A^/^t^/^ ast^oo. 

Remark 1.6. Compare with the Brownian case 

Px{TD{B)>t) r^C{x)e-^''^ ast^oo. 

Finally, let us point out the search for probabilistic connections with 
higher-order PDEs has gone on for a long time (at least since the 1960s). 
Biharmonic functions and related boundary value problems were studied 
by Vanderbei [47] and Helms [25, 26]. Concerning fourth-order operators, 
M§drecki [35, 36] explored connections with the 4-asymptotically stable 
motion and Nishioka [38, 39, 40, 41] considered the biharmonic pseudo- 
process — the one whose "transition density" is the fundamental solution of 
the parabolic biharmonic operator. Orsingher and his collaborators [1, 6, 7, 
37, 42, 43, 44] studied processes governed by signed measures related to fun- 
damental solutions of higher-order parabolic equations. Also, Hochberg and 
Orsingher [27] considered hyperbolic equations. Finally, Allouba [3] looked 
at connections with higher-order operators and Feynman-Kac-type formu- 
las. 

The article is organized as follows. The proof of Theorem 1.1 is given in 
Section 2. In Sections 3 and 4 the proofs of Theorems 1.3 and 1.5, respec- 
tively, are given. The last section is an appendix containing some technical 
results used in Sections 3 and 4. In Section 5 we prove Theorem 1.2. 

2. Proof of Theorem 1.1. Let / € C^(M"). As pointed out above, differ- 
entiation under the integral is not valid; hence to prove Theorem 1.1, we use 
distributions. Since u is bounded, it yields a distribution on O = (0, oo) x M"". 
Standard notation for this is u E P'(r2). We show u is a weak solution of 
(|A^ — -^)u{t,x) = — ^-^7=A/(x) on 0. Since the operator |A^ ~ ^ is hy- 

poelliptic and t~^/^A/(x) G C°°(r2), by the hypoelliptic regularity theorem, 
u S C°°{i}) and u is a classical solution to the partial differential equation. 
The initial condition is clearly satisfied and Theorem 1.1 will follow. 
To show u is a weak solution, it is enough to prove 

(2.1) / u{t,x) 

Jn 

for ip £ Co°(r2), the space of smooth functions with compact support in 
We describe the case n > 2, leaving n = 1 to the reader. Let e > 0. By 



-A^ + TT H 

8 dt 2^/2^ 



(p{t, x) dxdt = 
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Fubini's theorem, 



/ uA'^ifdxdt 



(2.2) 



Pn{y, \w - x\)f{w)pi(t,y)A'^ipdxdydwdt 



+ 2 



In Jo JB,_{wY 
= Ii + /2, say, 
where B^{w) = {x ■.\w — x\ < e} . Now 

^xPn{y, \w - x\)=Pn{y, \w - x\) 



Pn{y, \w - x\)f{w)pi{t,y)A'^(pdxdydwdt 



\w — x\ 



and, if denotes the outward normal derivative on dB^{wy, 

d e 
- — pn{y-, \w - x\) = -Pn{y,e) for X G dBs{w). 

OUa: y 

Hence by Green's identity, 

{A^ip)pn{y, \w - x\)dx 



(2.3) 



Pn{y,\w-X\) 



\w — x\ 



y^ 



Aipdx 



+ 



d e 
Pn{y,e)- — Aip- -pn{y,e)Aip 
OTix y 



da{x), 



ldB^{w) 

where (t{x) is surface measure on dBe{w). Another apphcation of Green's 
identity is not useful at this stage. Instead, observe 

§^P^(y^ 1^ ~ ^1) = \^wPn{y, \w - x\) 



\w — x\ 



Pniy,\w-X\). 



Hence, (2.3) becomes 

/ pn{y,\w-x\)A'^ipdx = 2 [ ■^[pn{y,\w - x\)]Aipdx 
JbawY JbawY oy 



+ 



dBeiw) 



Pn{y,e) 



— Aip--Ayj 
orix y 



da(x). 



Now we can write I2 from (2.2) as 

(2.4) /2 = /3 + /4 + /5, 
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where 

h = 2 



In Jo 
h = -2e 

h=4: 



n Jo 

oo 



f{w)pi{t,y)pniy,e 
f{w)pi{t,y)pn{y,e) 



d 



1 



Aif da{x) 



dy dw dt, 



y 



Aipda{x) 



dydwdt, 



f{w)pi{t,y) 



■^Pn{y,\w -x\ 
dy 



Aipdx dydwdt. 



In Jo JBe{wY 

By Fubini's theorem and integration by parts in the dy integral, 



1,^ = 4 



y 

f{w)pi{t,y)-pn{y, \w - x\)Aipdxdy dwdt. 
nJo JBe(wY t 



This has a form like I2, where now Aip in I2 is replaced by and pi(t,y) is 
t 



replaced by ^pi{t,y). Repetition of the derivation leading to (2.4) yields 



(2.5) 
where 



h = h + h + h, 



/6 = 4 
Ij = - 



nJo 
4e 



f{w)pi{t,y)-pn{y,s) 



dB^w) orix 



dip 



da{x) 



n Jo 

00 



f{w)pi{t,y)-pn{y,e) 



dBe{w) 



ipda{x) 



f{w)pi{t,y)^ip 
nJo JBe{wY t 



d_ 

dy 



Pn{y,\w-x\) 



dy dwdt, 
dydwdt, 
dx dy dw dt. 



By Fubini's theorem and integration by parts in the dy integral 

■1 y' 



nJO JBeiw) 



t t^ 



16 



nJO JBe(w) 



g-^m{t,y) 



pi{t,y)ippn{y, \w - x\)dx dydwdt 
ippn{y, \w — x\) dx dy dwdt 



using that ^pi{t,y) = — ^pi{t,y)[j — Now we integrate by parts in the 
dt integral to get 



(2.6) h 



-16 



°° f dip 

/ f{w)pi{t,y)—pn{y,\w - x\) dx dy dw dt. 
nJo Jb^(wY ot 



Since (p G Co°(r2), there is a bounded function h[t,w) with compact support 
in Q. such that 

d 

* '■ ' <h{t,w), xedBsiw) 



dn. 



-A'p>{t, x] 



and h is independent of e < 1. 
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Then for some Mi, M2, M3 > 0, 



|-^3|<C'/ / h{t,w)pi{t,y)pn{y,e)a{dBs{w))dydwdt 
Jn Jo 

jn-1 



<Ce" M / / pi{t,y)pn{y,e)dydwdt 

J Ml JBm,^{0)J0 

= Ce^-^ r" r^e-y'/' \ \ ,^ e-'"/'ydydt 

<Ce"-^ Z"*^' r ^ e-y'l'^^i^ \ e-^'^l'^y dvdt 

Jmi Jo V2^ (27ry)«/2 

Jo 



2 

The charge of variables u = |- transforms this to 



2?/ 

du. 



\l3\<Ce u''/^-^e^p(^-u- 



SMou 



If n = 2, 



/ u ^expf— M — — ^)(iti< / u ^expf — ^)dn+ / u^e^du 

Jo \ sNhv?) -Jo V sivwy A 



2 Je^/8M2 

after changing variables f = §1:^^^" in the first integral. Using this it is easy 
to show /3 — > as e — > 0. If n > 3, 

Then again, /a — > as e — > 0. A similar argument applies to Iq and Jy. In 
any event 

(2.7) hm Is = lim Iq = hm = 0. 

e^O e— ♦O e— ♦O 

By dominated convergence in (2.6), 

(2.8) lim Is = -8 [ u{t,x)^{t,x)dxdt. 
e-^O Jfi at 

Also, dominated convergence gives 

(2.9) lim h = 0. 
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/ uA"^ (f dx dt = lim( h + I2) 
Jn e^o' 



(2.10) 



lim lo 



lim(/3 + h + lQ + Ir + l8 

f d(p 
lim I4 — 8 u-?r- dx dt 



dt 



[using (2.4)-(2.5)] 
[by (2.7)-(2.8)]. 



It remains to analyze 1/^. Make the change of variables n = |^ in the dy 
integral to get 

2ei-" 



h 



(2.11) 



vr 



n/2 



I r\f{w)pikp)u-l^-^e-^ j A^da{ 
JnJo I \ 2u/ JdB^(w) 



dudwdt. 



Since f S Cq°{Q), there is a bounded function hi{t,w) with compact support 
in Q such that for e < 1, the dudwdt integrand (the portion square brackets) 
is bounded by 

Ct-^/\^/^-^e-''h{t,w)e''-\ 

where h and C are independent of e < 1. Then we can let e ^ in (2.11) 
and use dominated convergence to end up with 



lim Ia 



W2 jQJo 



1 



'2-Kt 



lim e 



1-n 



A{pda{x) 



du dw dt 



'2TTt 



27r"/2 



Lr(n/2) 



Aip{t, w) 



du dw dt 



-4 / f{w)—^A(p{t,w)dwdt 
JQ v2vrt 



1 



In v2vrt 
Hence, (2.10) reduces to 



Af{w)(p(t, w) dw dt. 



[ uA^ipdxdt = -A [ -^=Af{w)ip{t,w)dwdt-8 [ u^ 
Jn Jn V27rt Jn ot 

which is (2.1). 



dx dt, 
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3. Proof of Theorem 1.3. If D C M is an open set, write 
r^(x) = inf{t >Q:Xf + x^D}, 
and if / C M is an open interval, write 

r^l = r^{I)=mi{t>Q:YtiI}. 
By continuity of the paths, 

Px{TD>t) 

= P^{Zs£D ioT all s<t) 

= P{x + X+(0 V Ys) G D and X + X-(0 V (-F,)) G D for all s<t) 

(3.1) 

= P(r+(x) > V and (x) > V (-Y,) for all s<t) 
= P{-T^ (x) <Ys< T^{x) for all s<t) 
= P{ri{-T^{x),T+{x))>t). 

Once we show for n-dimensional Brownian motion B and a generalized cone 

(3.2) ^P.(rz)(i?)<t)~c(x)t-P(^)-i ast^oo, 

then Theorem 1.3 is an immediate consequence of (3.1) and the next theo- 
rem. 

Theorem 3.1. Let be a positive random variable with density f{t) 
such that for some positive C and p, f{t) ~ Ct^^~^ as t ^ oo. If and ^2 
are independent copies of ^, independent of the Brownian motion Y , then 
as t ^ 00, 

P{V{~^^,i;,) > t) ^ { t-Hnt, p=l, 

We give the proof of the asymptotic (3.2) in Lemma A. 3 of the Appendix. 

Proof of Theorem 3.1. We wiU abuse notation and allow Yq = x; 
the probability associated with this will be Px- For typographical simplicity, 
write 
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(3.3) Po(r?(_.,.) > t) = - E ^ exp 2(, + ,)2 V ■ 



The following distributions of t,) and r are well known: 

(2n^\fTi^ \ (2n + l)7rii 

(3.4) p^(r>t) = P,(r/(o,oo)>t) = 4= r'^\-^"l''dw 

\/2tt Jo 

([22], pages 340-342). By (3.4), 

(3.5) Px{V{o,i) >t)< P.ir > t) < ^ A 1. 

By Lemma A.l of the Appendix, choose M > so large that 

-^^^(^(O 1) > *) ~ (sin7rx)e~'^^*/^ for t>M, uniformly in x G (0, 1). 

(3.6) 

Then choose < | so small that 

(3.7) sinvrxssx, x£{0,6]. 
Choose K > so large that 

(3.8) f{u)^u-P-\ u>K. 
Finally, we consider t so large that 

(3.9) ''<'{i- 

By independence of y, , and ^2 > using scaling and translation invariance 
of Brownian motion, 

^(^(-?i,6)>0 

= ^o(??(-5i,6) >*) 

Po{V{-u,v) >t)f{u)f{v)dvdu 



JO 

00 roo / ^ 



, Pu/(u+v) V{o,i) > 7 — : — To f{u)f{v) dv du. 
Since Pxi^o,!) >t) = Pi-x{fl{o,i) >t),we have 

^ Pu/{u+v) ('7(0,1) > ^^l^y ) /(^)/(^) dv du 

Pv/iu+v) (V(0,1) > ^^l^y )fiu)fi'') dvdu 



00 fU 



JO 

00 roo 



10 Jv 

(reversing the order of integration) 



Pv/{u+v) (v{o,i) > J^^^^-^f{v)f{v)dudv 
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oo poo / -j- 

Pu/{u+v) ( ^(0,1) > 



/O Ju 

(relabeling). Thus, 



{u + v)' 



f{u)f{v)dvdu 



(3.10) 



2 / / Pu/{u+v) ( ??(0,1) > 



We break up the integral into several pieces: 



f{u)f{v)dvdu. 



Al = \{u,v)■.K<u<5^j^^,u + v>^|^ 



{u,v):Q<u<K,u + v>\l— \, 



^3 = I {u, v) : 5\j^ <u<v,u + v> 



{u,v) :u> 0,v > K,u <v,u + V < 



A^ = {{u,v):u> 0,u<v< K}. 
It turns out the integral on Ai is the "dominant" piece. 



The dominant piece: upper bound. On Ai, ^^^^ < ^^^^^[ = b and so > 



ItjM 



> ^^jj^u > K. Hence, by (3.5) and (3.8), 



t 



(3.11) 



< JJ^ ^^f{u)f{v)dvdu 



f{u)f{v)dvdu 



-1/2 



-1/2 



&Jt/M roo 



K 

5^/t/M 
K 
S^t/M 



't/M-u 



^P~^u~Pdvdu 



u du 



/ u-PiViydu (t large) 

JK 
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w ^ dw 



'K/Vi 

(changing variables u = \/tw). 



The dominant piece: lower bound. Notice on ^i,;^^^ < ^Y^*^.^ = 6 and 



/t/M 



< M. Hence, 



Pu/{u+v) [ ^(0,1) > 



f{u)f{v)dvdu 



> 



{u + vy , 

Pu/(u+v) (??{o,i) > M)f{u) duf{v) dv 



5^t/M poo 
K 



U 



t/M~u U + V 



-P-\-P-Uvdu [by (3.6)-(3.8)] 



S/VM roo 



W 



{z-w)-P-^w-P-Uzdw 



iK/Vt Ji/Vm z 
[changing variables z = (u + v)/ \/t, w = u/\/t] 



>t-P 



S/VM roo 
K/Vt Jl/y/M Z 



1 



dz 



w P ^ dz dw 
dw 



K/Vt 



z-P-^ dz 



l/VM 



rS/VM 

/ w~P dw for large t. 

JK/Vt 



Combined with (3.11), we end up with 



(3.12) 



u/{u+v) I V{0,1) > 



f{u)f{v)dvdu 



5/VM 
K/Vi 



w P dw. 



The remaining terms: upper hounds. On the set A2, v>\Jjj— K>K 
if t is large enough. Hence, by (3.5) and (3.8), 

// <r/" ^v-p~uvf{u)du 

J JA2 Jo JjlJU-u Vt 



K 



t-^l'^u 



u j{u) du 
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(3.13) 



« t-(p+i)/2 for t large. 
On the set A3, if t is large, then u > u > so by (3.5) and (3.8), 

<// 1- f{u)f{v)dvdu^ f f u-P-\-P-Uvdu 

A3 J J As J JA-i 

{Jt/M)/2 roa 

S^t/M Jy/t/M~u 

(3.14) + / / v-P~\-P-Uvdu 

Jiy/t/M)/2 Ju 

JSy/t/M p\2\ M J J{^/t/M)/2P 
w t-P + t-P 

^t-P. 

On A4, since t^jx^^t > M, by (3.6) and that ^i^^ is bounded for all x, 



A4 



Pu/{u+v) (^(0,1) > (^^^^2 ) /(^)/(^) 



AiU + V 

' -e-^"'/^^^+''^" f{u)v-P-^dvdu 



lo JK U + V 

-^''/^^^+-^\~p~\rP-Uvdu 



^tjM/2 rV^-u u _.2wo^,^..^2 



K Ju U + V 

[using (3.8)]. Changing variables z = w = u in the first integral and 



z = , w = in the second to get 



K fXlVU in 9, 9 

e-^ f{w)izVi-wrP-^Vidzdw 



J(w+K)/Vi Z\/t 



+ / / -e-^"/^'\wVtrP~\zVt-wVlrP~hdzdw 

JK/Vi J2w Z 

K 1-1/ VM ^ 9,^9^, . W \ -P-^ 



t-^P+^y^ / ^e-W/(^)^-P-i dzdw 

Jo J{w+K)/Vt Z \ zytj 
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+ t^P / -e"'' ^^"^ w Pz P ^{l dzdw. 



K/Vi J2w Z \ Z 

In the first term, since <w < K and < z, we have 

w w 1 
< rr < -: 



zVi~ w + K- 2' 

hence, (1 — < i'^) '^^^- the second term, since 2w < z, we have 



w 1 

7<2' 



hence, (1 - f^^^^ < i^yP-\ 
It follows that 



wz~P~^e~^ f{w)dzdw 

,w+K)/Vi 



+ t-P / z-P-^vj-Pe-'^/^'dzdw 

JK/Vi J2w 



(3.15) 



<t-iP+m / / wz-P-^e-^ f{w)dzdw 



JO 



rPi / z-f-^e"" / vj-Pdw 

/O / JiT/v^ 



/■1/(2VM) 



^-(p+l)/2^ // +t-P, 



using (3.12). 

Finally, on A^, since ^^^^^-^^ > M, by (3.6), 



As Jo JO 

(3.16) 

< t ^ for t large. 
Using (3.10) and (3.13)-(3.16), for t large, 

+ t-^P+^y^+t-P 



Ai 



< 
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By (3.12), for t large, 

/ r r \ f ^' 

// )^{rHnt, p=l, 

Hence, the desired conclusion follows. □ 

4. Exponential decay. In this section we prove Theorem 1.5. First we 
state a couple of lemmas. 

Lemma 4.1. Let ^ be a positive random variable such that for some 
c> 0, — logP(^ > t) ^ ct as t ^ oo. Then for independent copies and (,2 
of ^, — logP(^i + ^2 > ~ ct as t ^ oo. 

Proof. By independence, 

m + 6>t) = l- f P{i<t-y)dyP{i<y) 
Jo 

= 1- f\l-P{i>t-y)]dyP{i<y) 
Jo 

= P{C>t)+ f P{i>t-y)dyP{i<y) 
Jo 

>P{i>t). 

For an upper bound, note for any d < c, E[e^^ < oo. Then by the Markov- 
Chebyshev inequality and independence, 

Pi^i + ^2>t)<e-''E[e'^^'-^^'^] 
= e-''[E[e'i]f. 

Thus, 

et - 21og£;[e^«] < - logP(ei + 6 > i) < -log^(e > 
Divide by ct, let t ^ oo, then let | c to get the desired conclusion. □ 

We use Lemma 4.1 to derive asymptotics of the Laplace transform of 
(Ci +^2)"^- This is a simple application of de Bruijn's Tauberian theorem 
([9], page 254). We state a special case of it for the convenience of the reader. 
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DE Bruijn's Tauberian theorem (Special case). Let X be a positive 
random variable such that for some positive B and p, 

- \ogP{X <x)r^ Bx~P as x^O. 

Then 

-logEe"^^ ~ (p+ l)Bi/(p+i)p-p/(p+i)AP/(p+i) as X^oo. 

Lemma 4.2. With and ^2 as in Lemma 4.1, 

- logE{ exp(-A/(ei + 6)^)) ~ 3(c/2)2/3a1/3 as A ^ 00. 

Proof. By Lemma 4.2, 

- log p (^^^-^^ < = - log p(ei + 6 > 

— 1/2 

~ ex ' as X ^0. 
The claim follows from de Bruijn's theorem. □ 

We will also need the following application of de Bruijn's theorem. 

Lemma 4.3. Let X be a positive random variable with density f{u) = 
7u~^e~"/^ for u> 0. Here a > and 7 is chosen to give total mass 1. 
Then - logEe'^^ ~ 3a'^/^2-'^/^X^/^ as 00. 

Proof. We have, changing variables u = v~'^, 

P{X <x)= 7u-2e-"/^ du 
Jo 

POO 

= / 2jve-''"dv 

Hence, 

— logP{X<x)^ax~^^^ as X ^ 0. 
The desired conclusion follows from de Bruijn's theorem. □ 



The main result of this section is the next theorem. Before proving it, we 
show how it yields Theorem 1.5. 
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Theorem 4.4. Let > be a positive random variable whose density f 
satisfies — log/(t)~ct as t— >cx). //^i and ^2 o'^e independent copies of ^, 
independent of the Brownian motion Y , then 

- logP(7?(_5^,52) > t) ~ fc^/^t^/^ as t^ 00. 

In Lemma A. 4 of the Appendix we prove for a bounded open set D C 
with regular boundary, 

log(J^PArD{B)<t)^ r^-Xot ast^oo, 

where A/j is the principal eigenvalue of on with Dirichlet boundary 
condition. Then by (3.1) and Theorem 4.4, Theorem 1.5 follows. 

Proof of Theorem 4.4. First we give a lower bound. Let e > 0. Recall 
we use the notation f to mean for some constant ci > 0, g > cif . Choose 
M and 5 as in (3.6) and (3.7). By hypotheses, choose K >Q such that 

(4.1) e-<'^+e> < f^u) < e-^(i-^)" for u > K. 

For t large, let 

A=\^{u,v)■.K<u<5^|^^-^u<v<^|^^-u^. 

Then by (3.10), 

> ^ j Pu/(u+v) > (^y^l^y ^ fiu)fiv) dv du. 

On the set A, since ^ <\-, we have u > (^ — l)u > u> K and u + v > 
Thus, ^ < (5 and so by (3.6), (3.7) and (4.1), 



J{l~5)u/5 U + V 



(changing variables x = u + v,z = u) 



/ ^^-^Hl2x\-c(\+e)x ^^^^ 



IK Jz/5 X 

(reversing the order of integration) 
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Jk/s Jk X 



IK/S 

K/5 X 

> I _e-'^*/2^ 6-^(1+^)^(52x2 (fort large 

hx/s X 



> / ^^-.^t/2.^g-c(l+e)x ^2^2 _~__]dx 



^ 1 _,2.;o.2 
/2A7<5 X 

(changing variables u = x~'^ 

JM/t 

Thus we have for t large, 

(4.2) ^(r/(-a,6) > ^ / u-^e-<'+'y^e-- '"/^ 



Now 



(4.3) 



u 



U 



-2g-c(l+e)/v^^^<g-cit 



52 /4i^^ 



and changing variables u = v 2 

f-A//i 



(4.4) 



2g-c{l+e)/v^g-7r2tV2 



Jo 



9 r 

: z / ve 

/t/M 



By Lemma 4.3, with a = c(l + e) and A ■ 



2 ' 



(4.5) 



log / ^-2e-c{l+.)/v^g-^^t«/2 _ 3^2/3^2/3(1 ^ e)2/3il/3 

as t 
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By (4.3)-(4.5), we can rewrite (4.2) for large t as 

(4.6) PiV(-^u^,) >t)> exp(-|c2/3^2/3(i + ^)5/3^1/3^ 

Now we give an upper bound. By (3.10) and (3.6), 



^(^(-Ci,6)>*) = 2 



.+ 



U+V<y/t/M 

X ( Pu/{u+v) ( ^(0,1) > 



u+v>y/t/M 
t 



{u + vy 



f{u)f{v)dvdu 



(4.7) 



< 



u+v<Wt/M 



+ 



<E 



u+v>Jt/M 



J{u)f{v) dv du 



exp 



2(6 + 6)VJ^^«'^«=^V'/'''' 

Since > t) = Jt^ f{u) du, by our hypotheses on the density / of ^, 

log P{£, > t) ~ —ct as CO. 
By Lemmas 4.1 and 4.2, for t large (4.7) becomes 

^(^(-6.6) > ^ exp(-|c2/3vr2/3(l - ,)tV3) 



+ exp(-c(l -e)^f/M) 
<exp(-|c2/3vr2/3(l-e)ti/3). 

Combined with (4.6), this yields 

-|c2/37r2/3(l +e)5/3 < liminf i-i/3iQgp(^^_^^^^^^ > ^) 

< limsupt"^/^logP(r?(_^^_^2) > t) 

Let e ^ to get the desired conclusion. □ 

5. Proof of Theorem 1.2. In one dimension, with the boundary condi- 
tions = /(O) = /(I) = /'(O) = /'(I), the biharmonic operator is known 
to have a complete orthonormal set of eigenfunctions {ipn'n> 1} in L^(0, 1), 
with corresponding eigenvalues {A„:n > 1}. In fact, the eigenvalues are of 
the form A = a^, where a is a positive solution of 

cos(a) cosh(Q;) = 1, 
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and the corresponding eigenfunction is 

c(a) [(sin(a) +sinh(a))(sin(aa:;) — sinh(aa;)) 

(5.1) 

+ ( cos(q;) — cosli(a))( cos(aj;) — cosh.{ax))], 
with c(a) chosen to make the L'^{0, 1) norm one. Moreover, 



(2n+ l)7r' 



4 



{0.2} A„ ~ asn— >oo. 

These facts can be found in [19] on pages 113-116. 

Lemma 5.0. The eigenf unctions are uniformly bounded: 

sup |93„(rc)| < oo. 

n,x 

Proof. Simple computation reveals the constant c(a) in (5.1) satisfies 
c(a) ~ 2e~" as Q — > co. Expansion of (5.1) reveals the only trouble comes 
from the 

— sinh(a) sinh(ax) + cosh(a) cosh(aa;) 

terms. However, these combine to give cosh(a(l — x)) and clearly c{a) cosh(a(l — 
x)) is bounded for x G (0, 1). □ 

From (3.1), since P{t^ i)(^) > = PxiV{o,i) > we have 

/•oo roo 

(5.3) P^(r(o,i)(^) > t) Po{r]^^u,v)>t)fix,u)f{x,u)dudv, 

where f{x,u) = —-^Px{f]{o,i) > ^) is the density of f?(o,i)- This density has 
two expansions — one good for large u, the other good for small u: 

°° / C2?i+l)2vr2 \ 

(5.4) /(3;,n) =27r^(2n + l)exp --^ u sin((2n + l)7rx), 

n=o \ 2 J 

/(x,n) = (27r)-^/2^-3/2 

(5-5) 

X [(x + 2A;)e-("+2^')'/2" + (l-x + 2A:)e-(i-"+2'=)V2-] 

fc=— oo 

[see [22], page 342, and [10], page 172, (3.0.2)]. These expansions can be 
used to prove the following lemmas. We omit the details. 

Lemma 5.1. Forg{t,x) = Pa;(T(o,i)(-^) > t), the derivatives are bounded 
for x in compact subsets of (0, 1), t > and n < 4. 
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Lemma 5.2. The function g{t,x) = Px{t(^o^i){Z) > t) satisfies the bound- 



ary conditions jrS = x = 0, 1 for n<l 



Now we prove Theorem 1.2. To get a contradiction, assume 
5(t,x)=P,(r(o,i)(Z)>t) 

satisfies the equation 

dt^^^d^' 

Taking the Laplace transform in t and using Lemma 5.1 to interchange 
integration and differentiation, this becomes 

(5.6) Xg{X,x) -l = a-^g{X,x), 



where 



oo 



ff(A,x)= / e-^'git,x)dt. 



Since g is bounded, so is g; hence for each A > 0, we can expand with respect 
to the orthonormal basis {ipn-n'> 1} of L^(0, 1), 

oo 

g{X,x) = ^Cn{X)ipnix) in (0,1), 

n=l 

where 

Cn(A)= / g{X,x)ifnix)dx. 







To determine c„(A), note by (5.6) 

Acn(A) - / Lpn{y)dy = a 



A 



(Pn{x) dx. 



By Lemmas 5.1 and 5.2, -^g{X,x) =0 for x = 0,1 and n < 1. Since (fn 
satisfies the same boundary conditions, we can integrate by parts four times 
to get 

Ac„(A) - / ipn{y)dy = a [ g{X,x)(p'f> [x) dx 



-a / g{X,x)Xn'fn{x)dx 
Jo 

-aXnCn{X). 
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Thus, 



Consequently, 



Cn(A) 



Io'Pn{y)dy 
A + aXn 



g{X,x) = 2^ — X m L 0, 1 . 

By (5.2) and Lemma 5.0, the convergence is also pointwise. Then by the 
continuity theorem for Laplace transforms, for < a <b, 







(5.7) ^ / / ^n{y)dy 

For t large, say t>T 

oo r /.I 



LPn{x)e "■^"''ds= / g{s,x)ds 



E 

n=l 







fn{y)dy 



ipn{x)e 



-aXnt 



Vi{y)dy 
C{x)e-''^'' 



aXit 



^i(^)e-"^^ 



[by Lemma 5.0 and (5.2)] and 

exp(-i-iA(o('i)t'/')<5(t,x) 
(by Theorem 1.5). Hence for t>T, (5.7) yields 



{-cis^^^)ds< / e-^^i^ds. 



2t 



2t 

exp 
This implies 

Take the natural logarithm of both sides, divide by t^^^ and let t ^ oo to 
get 

-ci2^/3 < -oo, 

a contradiction. 



APPENDIX 
We collect some technical results used above. 

Lemma A. 1. Ast^co, 



4 2 

Px{V(o 1) > ~ — e"'^ sinvrx uniformly for x G (0, 1). 

TT 
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Proof. In light of the eigenfunction expansion 

A °° 1 
7r^j2n+l) 

([22], page 342), it suffices to show for some c> 0, independent of n, 

I sin(2n + l)7rx| , .9 , . 

' ^ ^ -<c(2n+l)2, j;E(0,l). 



sm vrx 
Since 

I sin(2n + l)7r(l — x)! | sin(2n + l)7rx| 



sin7r(l — x) sinvrx ' 

we can restrict attention to x G (0, There exists c> such that sinvrx > 
cx, X € (0, |] and so 



I sin(2?i + llvrxl (2n + l)7rx n , , /I 

<^ ^ — = -(2n+l), xgIo,- 



sm vrx cx c 



□ 



To prove (3.2), we collect some facts. If D C M" is a generalized cone, 
there is a complete set of orthonormal eigenfunctions rrij with corresponding 
eigenvalues < Ai < A2 < • • • satisfying 



Lgn-ivrij = —XjUij on D n S"' ^, 

mj = onS''~^ndD 

and for a{d9) surface measure on S*""^, 
[16]. For notational simplicity, set 



n ^2 



(A.l) ^^^V^^'^U"^ 
The confluent hyper geometric function (for 6 > 0) is 
^ , , , a z a(a + 1) z'^ 

Baiiuelos and Smits [5] have shown there is a series expansion 
(A.2) P.irniB)>t) = Y.B,{^-J^] i^^i 
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which converges uniformly for (x, t) € K x [T, oo), where K <^D is compact, 

r>o, 

(A.3) a,=j,-[^--l 
and 

(A.4) r((a,+n)/2) r (^)^(^^), 

^ ^ ^ r(aj+n/2) iDn5"-i ^ ^ ^ 

They also show for each x £ D, 

f X \ f IxP \ "^i/^ 

Hence, the p{D) described in the Introduction and appearing in the state- 
ment of Theorem 1.5 is ai/2. 

Lemma A. 2. Let k> be an integer. If a> k and 7 — a — 1 > 0, then 

< x"iFi(a,7, -x) < x^ ll^\l^''~^\ , x>0. 

r(a)r(7 — a) 

Proof. By formula 9.211.2 on page 1058 of [24], 

r(Q)r(7-Q) 



r(7) 



-iFi{a,'j,-x) 



x'" 



(-x)i~^ / eH''"\-x-ty-''-Ut 



7— cv— 1 



X 

'0 

X / \7-Q-l 

1 _ _ du 

\ xj 

X / \ 'y~Ci — 1 

e-"u"^i-V( 1- - ) du 
\ xJ 

< x^ r e-^'u"-^-^ du (since 7 - a - 1 > 0) 
Jo 

<x''r{a-k). 

It is clear from the nonnegativity of the integral in the third equality above 
that x"iFi(a,7,-2;) > 0. □ 

From formula 9.213 on page 1058 in [24], 

d (X 
(A.5) — iFi(a,7,x) = -iFi(q + 1,7 + 1,x) 

dx 7 
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and by Theorem 8, page 102, in [16], 

(A.6) sup <C(n)Aj"^^^/^ 

Also, the proof of Theorem 4.1 part (iv) in [20] carries through in our context 
to yield 



(A,7) 



< CXD. 



Lemma A. 3. The asymptotic (3.2) is valid. 

Proof. Let /(t) be the density of P:^{td{B) < t). Then 

d 



fit) 



dt 



PArD{B)>t). 



If we can show for fixed x the series obtained from (A. 2) by termwise dif- 
ferentiation with respect to t converges uniformly for t>T, then by (A. 5) 



\x\ 
~2t 



2\ aj/2 



n \x\ 



-i^i(->«i+2'-^ 



aj/2 
aj + n/2' 



n \x\ 



iFi l + ^,l + a,- + -,-— — 
i i I ^ o ' ^ J ^ o ' 2t 2t 



|2n 



Since iFi(a,7,z) — > 1 as z — > 0, the asymptotic (3.2) [with p{D) = ai/2] is 
immediate. Thus, we need only prove uniform convergence in t > T of the 
series in (A. 8) for x fixed. Choose so large that 



for a> N. 



T{a/2 - n - 4)(2a + n/2) ^ ^ ^ _ (.^3)/2 

r(a/2) 

Then choose J so large that 

aj> NM 2n for j > J. 

This is possible by (A.l) and (A. 3), since ^ oo as j — > oo by (A. 7). 
By (A.4) and (A.6), for t > T and j > J, 



BjlTlj 



^-1 



- ' J T{aj+n/2) 
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Hence, by Lemma A. 2 with k = n + 4, for j > J and t>T, the jth term in 
the series (A. 8) is bounded in absolute value by 



(n-i)/2 r((a-,- +n)/2) 
r(aj +n/2) 



2 \ n+4 



2 V 2t 



r(aj + n/2)r{aj /2 - re - 4) 
r(aj/2)r((aj+re)/2) 



+ 



1 ^\x\ 



2 \ n+4 



2 V 2t 



r(l + + re/2)r(l + aj/2 - re - 4) 



c(re)A} 



(n-l)/2 / |2; 



r(l + aj/2)r((oj +re)/2) 

2 \ n+4 



2t 



r(aj/2 - n - 4) {aj + n/2){aj/2 - re - 4)r(aj/2 - n - 4) 



r(a,-/2) 



+ 



(a,/2)r(a,/2) 



lirj r(a,/2) 



< 



2 



re 

2ao + - 



c(re)A 



(n-l)/2n2; 



2 \ n+4 



2t 



2"+6[aj(oj+re-2) 



-{n+3)/2 



by (A.9). 



Since aj{aj + n — 2) = Aj, this is the same as 



-c(re)2 



n+6 



\X\ 

~2t 



2 \ n+4 



X~'<c{n)T-'2 



lr)n+5 



2r 



2 \ n+4 



The uniform convergence for i > T of the series in (A. 8) follows from (A. 7). 
□ 

Lemma A. 4. For a hounded open set D C with regular boundary, 

\og(^j^P.^{TD{B) <t)^ ^ -Xnt ast^oo, 

where Xd is the principal eigenvalue of on D with Dirichlet boundary 
conditions. 

Proof. There is an eigenfunction expansion 

CO „ 

where < Ai < A2 < • • • and {^k} are the eigenvalues and eigenfunctions of 
on D with Dirichlet boundary conditions. Moreover, for some a > 0, for 
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each to > 0, 

(A.IO) |¥^fc(x)|2<at-"/V'=*° 
and 



oo 



E 



k=l 

([45], pages 121-127). Then it is easy to show for x £ D, 



oo „ 



and then that 
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